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Disorder effect on magneto-transport on the surface of a topological insulator 
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We study the magneto-transport properties on the disordered surfaee of a topological insulator attached to a ferro- 
magnet/ferromagnet junction. Since, in the surface Dirac Hamiltonian, an out-of-plane magnetization induces a mass 
gap, while an in-plane magnetization has a role of the effective vector potential, the mechanism of magneto-transport is 
different between these two cases. The former is similar to the conventional one in ferromagnetic metals, while the latter 
is due to the shift of Fermi circles in momentum space. Our numerical calculations show that the magnetoconductance 
in the in-plane configuration is robust against disorder compared to that in the out-of-plane configuration. 


Because of many possibilities of application to devices, a 
magnetoconductance effect has been searched for during the 
past several decades. In particular two-dimensional magneto¬ 
conductance effect is important to minimize devices. The lo¬ 
calization effect, however, is stronger in two dimension than 
that of in three dimension and thus the magnetoconductance 
is fragile against disorder.^ Indeed when time-reversal sym¬ 
metry is broken, all wave functions are localize in two dimen¬ 
sion but one exception at the critical point of the quantum Hall 
(QH) transition. 

Although extended wave functions are necessary for the 
magnetoconductance effect, the QH effect occurs in a strong 
magnetic field, and thus it is not practical to consider the mag¬ 
netoconductance in the QH regime. Recently, however, the 
QH effects without an external magnetic field have been real¬ 
ized on the surface of magnetically doped topological insula¬ 
tors (TIs).^’^ 

TIs are new quantum states of matters, which cannot be 
adiabatically connected to conventional insulators. A three- 
dimensional (3D) TI has a finite gap in the bulk but possesses 
gapless surface modes described by the two-dimensional (2D) 
massless Dirac Hamiltonian for simple cases.The surface 
states are dubbed as the helical surface, in which the spin 
quantization axis is perpendiculary locked to the momentum 
by spin-orbit coupling. 

The spin-momentum locking at the surface makes TIs 
promising for versatile device applications. At the interface 
of a ferromagnetic insulator and a topological insulator, a va¬ 
riety of unique magneto-transport phenomena has been theo¬ 
retically proposed^“^^ and experimentally examined. 

In this work we study the disorder effects on the mag¬ 
netoconductance of topological surface attached to a ferro- 
magnet/ferromagnet junction. The ferromagnetism on the sur¬ 
face is induced by the exchange interaction m • o", where 
a = (cTx, (Ty, (Tz) is the Pauli spin matrix of the surface elec¬ 
trons and m = (ntx, ruy, ruz) is the exchange field which has the 
direction of the magnetization and the magnitude of the ex¬ 
change splitting energy. The out-of-plane exchange field gen¬ 
erates a mass gap in the surface modes. When the Fermi level 
is located slightly above the bottom of the conduction band, 
by projecting into the conduction band, the surface states can 
be regarded as fully spin-polarized 2D electrons with conven¬ 
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tional parabolic dispersion. As in the conventional ferromag¬ 
netic metals, the conductance shows a change depending on 
whether the magnetizations of adjacent ferromagnets are in 
a parallel or an antiparallel alignment, namely magnetocon¬ 
ductance. On the other hand, the in-plane exchange field acts 
as an effective vector potential which shifts the Fermi cir¬ 
cles in momentum space. A misalignment of the Fermi cir¬ 
cles between two regions also gives rise to a magnetoconduc¬ 
tance.With the use of the transfer matrix method,we 
calculate magnetoconductance in the out-of-plane and the in¬ 
plane magnetization configurations, and compare the disorder 
dependence of them. Our result shows a difference between 
the two cases, that the in-plane magnetoconductance is rela¬ 
tively robust against disorder, compared with the out-of-plane 
magnetoconductance. 

We consider the surface of a TI which is attached to two fer¬ 
romagnetic insulators, FI and F2 (Fig. 1).^ The surface elec¬ 
tronic states in this system can be described by the 2D Dirac 
Hamiltonian 

H = V^PyCTjc - PxCTy) + X + U{x,y), (1) 

i=x,y,z 

where v is the velocity of the Dirac fermion, ruy and are 
the exchange fields and U(x,y) is the disorder potential. We 
note the exchage field in FI as mi = (mix, miy, m\z) and F2 
as m 2 = (m 2 x, m 2 y, m 2 z)- We assume the TFs surface which is 
attached to two ferromagnetic insulators as 2D sheet of length 


nil m2 



Fig. 1. Schematic picture of F1/F2 junction on the TI surface.^ The ferro¬ 
magnetism on the surface of a TI is induced due to proximity effect by the 
ferromagnetic insulators deposited on the surface. The current flows on the 
surface of a TI. 
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2L, where the length of each ferromagnetic insulator is L, 
along longitudinal x direction and width W along transversal 
y direction. We take the aspect ratio W/L = 4, 8,12. 

We utilize the transfer matrix method^^“^^ for the 2D Dirac 
equation if'P = where 'F(x,y) is the two-component 
(spinor) wave function. Multiplication of both sides by icTy 
gives 

fivd^ = {vcTzPy + i (U(x, y) - E)cry rUxCTz + irriy - rrizCTj^ 'P. 

( 2 ) 

Here we discretize variable v at the N points xi,X 2 , ...,xjv. The 
transfer operator At, defined by ^ 2 l = At'Fo, is given by the 
operator product 


V = 0, the sum of incoming and refiected waves given by 

k' 

while the sum of transmitted waves at v = 2L is given by 

= (12) 

k' 

The right moving component in mode k is 0^ and left moving 
component is 0^. Starting from a mode incident from right, 
we can similarly obtain the refiection and transmission matri¬ 
ces r' and t\ which give, with r and t, the unitary scattering 
matrix. 




Km ,Xi = exp [(1 /n) apyCr^ . (4) 

The operator V gives the decay of evanescent waves between 
two scattering events, described by operator 7C„, 





(5) 

where 


- exp ^(ia/fiv) (Un - E)cry ^, 

(6) 

and 

Kn 

= Qxp[(a/nv)mxcrz], 

(7) 



= exp |^(/a/liv)myj, 

(8) 


Kn 

= exp [(-a/fiv) rUzCrx] , 

(9) 

where a is the lattice constant. 


To calculate the transfer matrix, we represent the operators 

in the basis 




1+ 

II 

1 Ink 

^/w w 

(10) 

The spinors |+) = 

2 eigenvectors of 

-CTy. By 


truncating the transverse momenta qk at \k\ = M, the dimen¬ 
sion of the transfer matrix becomes finite. The disorder po¬ 
tential U(x, y) = Zn,m - Xn)S{y - y^) is given by a col¬ 

lection of isolated impurities distributed uniformly over the 
scattering region 0 < v < 2L, 0 < y < W. The strengths 
jnm of the scatterers are uniform in the interval [- 70 . To]- 
The disorder strength is quantified by the correlator ^0 = 
f dr'(U(r)U(r')) which evaluates to Kq = ^y^Cl/hva)^, 
independent of the correlation lengths. The disorder strength 
can be related to the mean free path ^ in Boltzmann trans¬ 
port limit by ^ = TivIKqE. The average conductance (G) is 
obtaibed by sampling some 200-2000 disorder realizations of 
the impurity potential. We take ILja large enough so that the 
calculation result does not depend on the orders of operator 
products {ILja > 20). The momentum cutoff M is also large 
enough in this calculation (M > 20). 

To formulate the scattering problem, we consider a 
scattering state that has unit incident current from the left 
(v = 0) in mode k and zero incident current from the right 
(x = 2L). The quantum number k labels transverse modes. At 



As a consequence of unitarity, the matrix product and 
have the same eigenvalue called transmission eigenval¬ 
ues. The conductance G follows from transmission eigenval¬ 
ues via the Landauer formula G = jTr 

The information contained in the scattering matrix S can 
equivalently be represented by transfer matrix At.^^’ While 
the scattering matrix relates outgoing waves to incoming 
waves, the transfer matrix relates waves at the right to wave 
at the left, 

xj/right _ 


We separate the spinor degree of freedom of At into four 
blocks 


At = 


/At^+ 

\At-+ 



(15) 


As one can verify by substitution into Eq. (14), and compar¬ 
ison Eq. (11) and Eq. (12), the submatrices At^^' (s, s' = +) 
are related to the transmission and refiection matrices by 

r = - M-^, r' = M"' (16) 

t = (A)-)"' M-^, t' = (M"")"' (17) 


The repeated multiplication of transfer matrices is unstable 
because it produces both exponentially growing and exponen¬ 
tially decaying eigenvalues, and the limited numerical accu¬ 
racy prevents one from retaining both sets of eigenvalues. We 
resolve this obstacle by converting the transfer matrix into a 
unitary matrix, which has only eigenvalues of unit absolute 
value. 

This completes the description of our numerical method. 
We now turn to the result. The inset in Eig. 2 shows the dis¬ 
order strength dependence of conductance for £" = 0 and 
E = {)21ftv I a in the absence of exchange field. The disorder 
dependence qualitatively changes around ~ 0.3, indicating 
the crossover between the ballistic and the diffusive regimes. 
Indeed the mean free path is comparable to the system size 
around this point. The main panel shows the magnetoconduc¬ 
tance in the out-of-plane exchange field configuration. Here 
we fix m 2 z = 0.2flvla and vary from -m 2 z to +m 2 ^. The 
Eermi energy E is also fixed at £" = 0.27fiv/a = 1.35m2^. The 
normalized magnetoconductance is plotted as function of 
for various disorder strength. 

In the clean limit, Kq = 0, the magnetoconductance is pos¬ 
itive for the parallel configuration (mizlm 2 z > 0), while it is 
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Eig. 2. The out-of-plane exchange field dependence of normalized conduc¬ 
tance for various disorder strength. Here mi = (0,0, mi^), m 2 = (0,0, m 2 z) 
and m\z is varied from -m 2 z to +m 2 z with fixing m 2 z = O.lfivla. The aspect 
ratio and the energy are fixed at W/L = 4 and E = 1.35m2z- The inset shows 
the disorder strength dependence of the conductance for E = 0 (pink) and 
E = O.llfivla (blue). 


negative for the antiparallel configuration (miz/m 2 z < 0). We 
see good agreement with analytical results of the out-of-plane 
field dependence in the continuum model.^ 

To understand this magnetoconductance behavior, we con¬ 
sider the case where |m^| » vtikp, kp being Fermi wave num¬ 
ber, and focus on the positive energy band on the surface of 
a TI. In this case, the spin degeneracy is lifted and the spin 
direction is (-fivky, fivkx, rriz). The energy dispersion approx¬ 
imated 2 iS E = ^(fivk)^ -f mj ~ -f |m^|, regarded as 

fully spin polarized Schrodinger electrons. For this reason, 
the mechanism of this magnetoconductance induced by out- 
of-plane exchange field corresponds to that in a conventional 
ferromagnetic metal. We note the Landau level gaps are negli¬ 
gibly smaller than the energy gap induced by the out-of-plane 
exchange field. 

Next, we consider the infiuence of disorder on the surface 
of a TI. Figure 2 shows that the dependence of the conduc¬ 

tance changes abruptly when weak disorder is introduced. The 
conductance {G(miz)) does not become maximum but mini¬ 
mum at mi^ = 4-m2z (parallel configuration), in contrast to the 
clean limit case. In the disordered case the conductance takes 
it maximum at ~ 0. At strong disorder, .^0 = 1-0, 
dependence of the conductance is nearly symmetric around 
miz = 0, indicating that the parallel/antiparallel configuration 
does not matter. 

When the exchange field is applied in the v direction, 
kyCTx term in the original Dirac Hamiltonian is replaced by 
(ky -f mxlfiv)(Tx indicating that the Fermi circle is shifted by 
-nix/fiv in momentum space, while the size of Fermi circle re¬ 
mains unchanged. When m\x ^ m 2 x, the position of the Fermi 
circles in FI and F2 are different. This misalignment of the 
Fermi circles causes a change of conductance (magnetocon¬ 
ductance) with qualitatively different mechanism from con¬ 
ventional one. Namely, as \m 2 x - m\x\ increases, the overlap 
region of the Fermi circles between FI and F2 is reduced as 



'mixlm2x 

Eig. 3. The in-plane exchange field dependence of normalized conduc¬ 
tance for various disorder strength, mi = (mi;c,0,0), m 2 = (m 2 x, 0 ,0) and 
mix is varied from -m 2 x to +m 2 x with fixing m 2 x = 0.2fivla. The aspect ratio 
and the energy are fixed, W/L = 4 and E = 1.35m2x- 


depicted in the inset in Fig. 3, and thus the number of evanes¬ 
cent modes increases, therefore the conductance decreases. 
The in-plane exchange field induced magnetoconductane is 
characteristic to the spin-momentum locking and qualitatively 
different from the conventional magnetoconductance. 

The normalized magnetoconductance in the in-plane con¬ 
figuration, mi = (mix, 0,0) and m 2 = (m 2 x, 0,0), is shown in 
Fig. 3. We fix the exchange field in F2 at m 2 x(ci/fiv) = 0.2, 
and vary the exchange field from -m 2 x to -\-m 2 x- In the 
clean limit (Kq = 0), the conductance increases for the par¬ 
allel configuration (m\xlm 2 x > 0), and decreases for the an¬ 
tiparallel configuration (m\xlm 2 x < 0). As mentioned above, 
the conductance is influenced by the relative positions of the 
Fermi circles. The conductance takes maximum when there 
is no misalignment of Fermi circles between FI and F2, 
and decreases with introducing the misalignment. Again we 
see good agreement with analytical result in the continuum 
model.^ 

As Fig. 3 shows, the dependence of magnetoconductance 
on the in-plane exchange field becomes gradually weak as the 
disorder strength increases, in contrast to the case of out-of- 
plane dependence. Even in the presence of disorder, the mag¬ 
netoconductance remains positive for the parallel configura¬ 
tion (mix/m 2 x > 0), while negative for the antiparallel config¬ 
uration (mixlm 2 x < 0). 

To compare disorder dependence of the magenetoconduc- 
tance with the out-of-plane exchange field and that with the 
in-plane, we plot the conductance differences {G(miz/x = 
m 2 z/x, Ko)} - {G(mizix = -m 2 z/x, ^o)> as a function of disor¬ 
der strength Kq in Fig. 4. The conductance difference with the 
out-of-plane field abruptly decreases as the disorder strength, 
.^ 0 , increases even in the ballistic regime and vanishes in the 
diffusive regime. On the other hand, the conductance differ¬ 
ence with the in-plane field decays more slowly and remains 
finite even in the diffusive regime. These behaviors are found 
at all the aspect ratios we examined, W/L = 4, 8,12. These 
results clearly indicate that the magnetoconductance with the 
in-plane field is robust against disorder while that with the 
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Eig. 4. The disorder dependence of the normalized conductance difference 
between parallel and antiparallel exchange field conhgurations for out-of- 
plane (red) and in-plane (blue) case. 


out-of-plane field is fragile. Since the former is characteristic 
to the surface states of a topological insulator (strong spin- 
orbit coupling), it is a great advantage of topological insulator 
based devices. 

In conclusion, we have studied the disorder effect on the 
magnetoconductance of the ferromagnet/ferromagnet junc¬ 
tion on the surface of a TI. With the use of the transfer ma¬ 
trix method, we calculated the magnetoconducrance in both 
the out-of-plane and the in-plane exchange field configura¬ 
tions. In the out-of-plane field, the Dirac electrons are re¬ 
garded as fully spin-polarized Schrodinger electrons when the 
Fermi level is located slightly above the bottom of the conduc¬ 
tion band. The mechanism of the magnetoconductance in this 
regime corresponds to that in a conventional ferromagnetic 
metal. On the other hand, the in-plane field induced magne¬ 
toconducrance is characteristic to the surface of a TI. These 
two cases show different disorder dependence. These results 
are consistent with the fact that all wave functions are local¬ 
ized in the presence of a mass gap, while in the in-plane fields 


massless Dirac fermion systems belong to the critical point 
of the quantum Hall transition, and thus wave functions are 
extended.^’ ^ Since the latter is robust against disorder, it is an 
advantage of TI based devices. 
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